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During development, multipotent progenitor cells 
establish tissue-specific programs of gene expres- 
sion. In this paper, we show that p63 transcription 
factor, a master regulator of epidermal morphogenesis, 
executes its function in part by directly regulating expres- 
sion of the genome organizer Satbl in progenitor cells. 
p63 binds to a proximal regulatory region of the Satbl 
gene, and p63 ablation results in marked reduction in the 
Satbl expression levels in the epidermis. Satbl ^ mice 
show impaired epidermal morphology. In Satbl -null epi- 
dermis, chromatin architecture of the epidermal differen- 
tiation complex locus containing genes associated with 



epidermal differentiation is altered primarily at its central 
domain, where Satbl binding was confirmed by chroma- 
tin immunoprecipitation-on-chip analysis. Furthermore, 
genes within this domain fail to be properly activated 
upon terminal differentiation. Satbl expression in p63 +/ ~ 
skin explants treated with p63 small interfering ribo- 
nucleic acid partially restored the epidermal phenotype 
of p63-deficient mice. These data provide a novel mecha- 
nism by which Satbl , a direct downstream target of p63, 
contributes in epidermal morphogenesis via establishing 
tissue-specific chromatin organization and gene expres- 
sion in epidermal progenitor cells. 



Introduction 

During the execution of developmental programs, the compo- 
nents of the chromatin-remodeling and transcription machinery 
are expressed in a highly organized manner, resulting in the 
establishment of tissue-specific patterns of gene activation and 
silencing (Hemberger et al., 2009). However, the mechanisms 
underlying these series of events to achieve terminal differenti- 
ation of specific cell types are still unclear. 

The program of epidermal differentiation and barrier for- 
mation in mice begins at about embryonic day 14.5 (E14.5) 
followed by the appearance of the spinous and granular layers 
at El 6.5 and establishment of a functional epidermal barrier 
at E18.5 (Fuchs, 2007; Koster and Roop, 2007; Blanpain and 
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Fuchs, 2009). This program is tightly controlled by several 
transcription regulators, including the p63 transcription fac- 
tor (Koster and Roop, 2007; Truong and Khavari, 2007; Crum 
and McKeon, 2010), as well as by the DNA- and chromatin- 
remodeling factors (DNA methyltransferase DNMT1, histone 
demethylase JMJD3, histone deacetylases 1/2, ATP-dependent 
chromatin-modifying enzymes Brgl and Mi-2p, and the poly- 
comb group protein Ezh2; Indra et al., 2005; Kashiwagi et al., 
2007; Sen et al, 2008, 2010; Ezhkova et al., 2009; LeBoeuf 
et al., 2010). It remains to be explored how transcription factors 
and chromatin-remodeling enzymes are coordinately regulated 
to establish tissue-specific gene expression programs during 
terminal differentiation of epidermal progenitor cells. 
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The p63 transcription factor serves as a master regulator 
of epidermal development, and p63 knockout (— /— ) mice fail 
to form stratified epithelium and to express several epidermis- 
specific genes (Mills et al., 1999; Yang et al., 1999). In kera- 
tinocytes, p63 regulates a large number of genes that encode 
distinct adhesion/signaling molecules, transcription factors, 
and cell cycle-associated proteins as well as tissue-specific 
proteins, such as keratins, involucrin, and loricrin (Vigano and 
Mantovani, 2007). 

Eukaryotic genomes appear to be organized for gene regu- 
lation with a tendency of linear clustering of coexpressed genes, 
often related in function (Kosak and Groudine, 2004). In the 
mammalian epidermis, differentiation of the multipotent pro- 
genitor cells in the basal epidermal layer into keratinocytes of 
the suprabasal layer is accompanied by coordinated activation 
of the specific sets of genes; e.g., a cluster of such genes is found 
in the epidermal differentiation complex (EDC) located in the 
gene-rich region of mouse chromosome 3 (Fig. SI). Genes in 
the EDC encode the components of the cornified cell envelope 
(involucrin, loricrin, filaggrin, small proline-rich and late- 
cornified envelope proteins, etc.) essential for epidermal barrier 
function (Marshall et al., 2001; Martin et al, 2004; Bazzi et al., 
2007; Brown et al., 2007). 

It is well accepted now that higher-order chromatin organi- 
zation and spatial arrangement of genes within the nuclear space, 
as well as nuclear compartmentalization of chromatin-remodeling 
complexes and transcription machinery, all play an important role 
in controlling gene expression (Fraser and Bickmore, 2007; 
Lanctot et al., 2007; Takizawa et al., 2008; Rando and Chang, 
2009; Zhao et al, 2009; Joffe et al, 2010; Naumova and Dekker, 
2010; Schoenfelder et al., 2010). It has been previously shown 
that the differentiation process in cultured epidermal keratinocytes 
is also associated with changes in intranuclear positioning of the 
EDC and distinct chromosomes (Elder and Zhao, 2002; Marella 
et al., 2009). siRNA knockdown of p63 in keratinocytes results in 
marked alterations of gene expression in the EDC (Truong et al, 
2006). However, it remains to be determined whether the p63 pro- 
tein itself directly regulates some of these events during epidermal 
differentiation or whether the proteins encoded by p63 down- 
stream target genes play a role in some of these processes. 

In several other tissue-specific gene loci (e.g., T H 2 cyto- 
kine locus and p-globin locus), higher-order chromatin remod- 
eling for establishing specific 3D conformations is regulated 
by the specialized adenine and thymine-rich binding protein 
Satbl, which functions as the genome organizer (Dickinson 
et al., 1992; Alvarez et al., 2000; Cai et al., 2003, 2006; 
Kumar et al., 2007; Han et al., 2008; Wang et al., 2009). Satbl 
targets chromatin-remodeling enzymes and transcription fac- 
tors to specific genomic regions, establishes region-specific 
epigenomic modification status, and plays a fundamental role 
in the execution of tissue-specific gene expression programs 
(Yasui et al., 2002; Cai et al., 2006; Pavan Kumar et al., 2006). 
However, the upstream mechanisms controlling Satbl expres- 
sion in distinct cell types are unknown. 

Here, we show that during epidermal morphogenesis, p63 
transcription factor directly regulates expression of the genome 
organizer Satbl. Satbl remodels chromatin architecture at the 



tissue-specific EDC gene locus, as indicated by compression of 
the central domain of this locus where terminal differentiation- 
associated genes are clustered. Upon Satbl ablation, these genes 
cannot be properly activated, and the epidermal morphology 
becomes impaired, thus demonstrating an essential role for Satbl 
as a part of the p63-dependent developmental program, through 
higher-order chromatin remodeling and gene regulation in epi- 
dermal progenitor cells. 

Results 

Expression of genes encoding chromatin- 
remodeling factors is markedly altered in 
the p63-null epidermis 

Because p63 transcription factor is a master regulator for epider- 
mal development, we hypothesized that it may regulate expres- 
sion of a multitude of genes encoding chromatin-remodeling 
factors to establish specific chromatin architecture in epidermal 
progenitor cells. Global gene expression profiles were com- 
pared between the epidermis isolated by laser capture microdis- 
section (LCM) from E16.5 p63~'~ and wild-type (WT) mice, 
as previously described (Fig. 1 A; Mammucari et al., 2005; 
Sharov et al., 2006). This approach allowed us to compare gene 
expression profiles between the epidermal cells of p63~'~ and 
WT mice in situ and to exclude alterations in gene expression 
induced by the proteolysis-associated cell separation proce- 
dures (Zhou et al., 2006). Global microarray analysis showed 
> 1.8-fold changes in expression of 2,522 up-regulated and 998 
down-regulated genes between the epidermis of p63 _/ ~ and WT 
mice (Fig. 1 B and Tables SI and S2). 

Among the distinct groups of genes, the expression of 
which was either down- or up-regulated in p63~'~ mice versus 
WT mice, gene ontology analysis revealed a significant enrich- 
ment of the genes (P < 0.05) involved in the control of nuclear 
and chromatin assembly and remodeling (Fig. 1 B). These genes 
constituted ^4 1 and 5% of the total number of down-regulated or 
up-regulated genes, respectively. Interestingly, in the epidermis 
of p63 _/ ~ mice versus WT controls, the expression of the genes 
encoding several key regulators of the higher-order chromatin 
structure and ATP-dependent chromatin remodeling (Satbl, 
Mi-2a, and Mi-2p) was down-regulated, whereas the expression of 
some other ATP-dependent chromatin-remodeling factors (Chdl 
and Brm) was up-regulated. The top list of most affected genes is 
shown in Fig. 1 B (see the entire list in Tables SI and S2). The 
relative expression levels of some of the most down-regulated 
genes between the WT and p63~ /_ mice were further validated by 
RT-PCR (Fig. 1 C). These data suggest that the genes encoding 
chromatin-remodeling factors represent a novel class of the p63 
targets, raising a possibility that some of these factors may be in 
charge of the remodeling of the higher-order and/or local chro- 
matin structures and tissue-specific gene expression during dif- 
ferentiation of the epidermal progenitor cells. 

p63 controls the expression of Satbl in 
epidermal progenitor cells 

We studied whether p63 controls tissue-specific chromatin organi- 
zation and gene expression during development of the epidermis 
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Figure 1 . Changes in global transcription profile and 
expression of the genes encoding chromatin-remodeling 
factors in the epidermis of p63~'~ mice. The skin of 
El 6.5 p63~/~ and WT mice was processed for LCM 
to isolate fragments of the epidermis. RNA was iso- 
lated from captured tissue, amplified, and processed 
for microarray and RT-PCR analyses. (A] Cryosection 
of the skin after LCM including a fragment of the cap- 
tured epidermis. Bar, 25 pm. (B) Agilent microarray 
analysis of the laser-captured epidermis of p63~^~ 
and WT mice. Diagrams showing the ontology of the 
down- and up-regulated genes in p63-null versus WT 
epidermis. Selected genes involved in the control of 
nuclear structure and chromatin remodeling, the ex- 
pression of which was changed in the epidermis of 
p63~/~ versus WT mice, are listed (a full list of the 
genes is shown in Tables SI and S2). (C) Real-time 
PCR for Serb J, Mi-2a, Mi-2/3, Terf2, and Prmt4 in the 
El 6.5 epidermis of p63~^~ mice normalized to the 
corresponding levels in age-matched WT mice. Error 
bars represent SEM. 
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through regulation of some of the chromatin-remodeling factors. 
For this purpose, we focused on the genome organizer Satbl as 
a potential p63 target because of its known function in regulating 
large-scale chromatin remodeling in several cell type-specific 
gene loci (e.g., T H 2 cytokine locus and p-globin locus; Cai et al., 
2003; Wang et al., 2009) as well as cancer-related gene loci in 
aggressive breast cancer (Han et al., 2008). We found that Satbl 
was coexpressed with p63 in basal cells of the epidermis of em- 
bryonic (E16.5) and postnatal skin (Fig. 2 A and not depicted). 
Skin epithelium of p63~'~ mice expressing Keratin 18 as a 



marker of undifferentiated epithelial cells showed strongly re- 
duced Satbl protein levels (Fig. 2, B and C), which was consis- 
tent with greatly reduced Satbl transcript expression in the 
epidermis of p63 _/ ~ mice versus WT controls (Fig. 1 C). 

In addition, the expression of the Satbl protein was 
strongly reduced in the E13.5 skin explants of heterozygous 
p63 knockout (+/— ) mice treated with p63 siRNA when com- 
pared with age-matched skin samples treated by the control 
siRNA (Figs. 2 D and S2 [A and B]). The expression of Loricrin, 
a marker of differentiated keratinocytes, was also decreased in 
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Figure 2. p63 controls the expression of Satbl in epidermal progenitor cells. The skin of p63-deficient and WT embryos was processed for immuno- 
fluorescent analyses as well as for modulation of gene expression in organ culture experiments with p63 siRNA. (A] Double immunofluorescence show- 
ing coexpression of p63 and Satbl in the nuclei of basal epidermal cells in E16.5 WT mice. (B and C) Immunostaining for Satbl (B) and K18 (C) in 
the El 6.5 skin of WT and p63~^~ mice. Marked decrease of Satbl expression in the epidermis (EP) of p63~^~ mice is magnified in the dashed boxes. 
(D) Immunostaining for p63, Satbl , and Loricrin in the El 3.5 skin samples of heterozygous p63 knockout (+/— ] mice cultured in the presence of p63 siRNA 
or control siRNA for 48 h. In epidermis, the expression of p63, Satbl , and Loricrin proteins in the p63 siRNA-treated samples was decreased compared 
with the controls. (B and D) The epidermal-dermal junction is shown by dotted lines. Bars, 25 pm. 



the epidermis of embryonic skin explants after treatment with 
p63 siRNA compared with controls (Fig. 2 D). 

Furthermore, chromatin immunoprecipitation (ChlP) assays 
using primary epidermal keratinocytes isolated from E16.5 em- 
bryos or newborn mice revealed binding of the p63 to the con- 
sensus sequence identified in the upstream region at —1662 
to —1738 bp from the Satbl transcription start site (Fig. 3, 
A and B). p63 binding was also seen at the consensus sequence 
in the Cldnl promoter region, an established p63 target gene 
(Fig. 3 B; Lopardo et al., 2008). However, no binding was seen 
at the predicted p63 site located in the first intron of the Satbl 
gene, serving as a negative control for this assay (Fig. 3 B). 

To assess whether p63 controls the activity of the Satbl 
promoter, HaCaT keratinocytes were cotransfected with a lu- 
ciferase reporter plasmid containing a mouse Satbl promoter 
(pSatblluc; unpublished data) and either the ANp63 expression 
vector (pANp63) or a control vector (pcDNA3). Transfection of 
cells with pANp63, but not with pcDNA3, resulted in an increase 
of the Satbl promoter-driven luciferase activity by threefold 
compared with control (Fig. 3 C). These data strongly suggest 
that Satbl is indeed a direct target gene of p63 in keratinocytes. 



pG3 or Satbl ablation results in altered 
conformation of the chromatin domain 
containing the EDC locus 

We further investigated whether there is any evidence for 
changes as a result of either p63 or Satbl ablation in the higher- 
order chromatin organization in epidermal progenitor cells dur- 
ing development. Using 3D FISH, we analyzed the conformation 
of the 5 Mbp chromatin domain of mouse chromosome 3 that 
contains the tissue-specific EDC gene locus in the epidermis of 
E16.5 p63 - ' - , Satbl and corresponding WT embryos. The 
EDC occupies ~3.1 Mbp in this domain located in the gene- 
rich region of mouse chromosome 3 and is flanked by several 
genes that are expressed in the epidermis, including Rps27 
and Gabpb2, encoding ribosomal protein S27 and guanine and 
adenine-binding protein subunit 2, respectively (Fig. 4 A). 

Despite the fact that the mean radius of the nuclei of basal 
epidermal cells was quite similar in p63 -/- , Satbl and cor- 
responding WT embryos, the distances between the Rps27 and 
Lor genes in the E16.5 epidermis of both p63 -/- and Satbl - ' - 
mice (without or after normalization to the corresponding nu- 
clear radius) were significantly increased (P < 0.05) versus the 
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corresponding controls, whereas distances between the Lor and 
Gabpbl genes remained unchanged (Figs. 4 [B-D] and S2 C). 
Distribution of the distances between the Rps27 and Lor in 
individual nuclei showed a lack of peaks in the corresponding his- 
tograms, thus suggesting similarities in the positioning of both 
alleles of each gene relative to other genes in the keratinocyte 
nuclei (Fig. S2 D). However, there was no significant difference 
observed in the distances between any of those genes in dermal 
cells of the p63 -/- or Satbl - ' - mice and WT controls, indicat- 
ing that the effects of p63 and Satbl on 3D chromatin structure 
in keratinocytes are tissue specific (Fig. 4 D). 

Next, we examined whether there is any correlation of the 
changes in chromatin organization and gene expression at the 
EDC locus for E16.5 epidermis of p63 - ' - and Satbl - ' - em- 
bryos. The EDC contains several genes, including those acti- 
vated during terminal keratinocyte differentiation (Fig. 4 A). 
Using quantitative RT-PCR, we observed that genes within the 
EDC and its flanking regions (Rps27 and Gabpbl) showed quite 
similar changes in expression in the epidermis of El 6.5 p63 - ' - 
and Satbl - ' - embryos versus the WT controls (Fig. 4 E). 

In both p63 - ' - and Satbl ~'~ mice, the expression of the 
genes located in the central EDC domain and involved in the 
epidermal barrier formation, such as Lor, Inv, Lcelb, and Lce3c, 
was decreased, whereas the expression of the S100a9 located at 
the EDC 5 ' end increased versus the corresponding WT controls 
(Fig. 4 E). In Safbl -/- mice, similar changes in gene expression 
in the EDC were observed in both dorsal skin and foot pads 
(Figs. 4 E and S2 E). The higher levels of reduction in expres- 
sion for some of the genes tested in the epidermis of p63 - ' - 
mice compared with those of the Satbl ~'~ mice at El 6.5 suggest 
the requirement of additional factors downstream of p63 for 
regulation of their expression. Expression of Fil, located at the 
EDC 3' end, is solely dependent on p63 and not on Satbl, indi- 
cating that Satbl is insufficient to replace global gene regulation 
by p63 (Fig. 4 E). Nevertheless, a similarity in altered expres- 
sion patterns of the EDC genes seen in the p63 - ' - and Satbl ~'~ 
mice, paralleled by similar changes in the conformation of the 
5 Mbp chromatin domain containing EDC, suggests a role for 
Satbl in mediating, at least in part, the p63 functions in the 
developing epidermis. 

Satbl ' mice have impaired epidermal 
morphology and altered gene expression in 
keratinocyte -specific genomic loci 

If Satbl is a critical downstream target gene of p63 in epider- 
mal progenitor cells, it is expected that the skin of Satbl - ' - 
mice will exhibit altered epidermal morphogenesis. In fact, 
Satbl - ' - mice at postnatal day 0.5 (P0.5) show a significant 
decrease (P < 0.01) of the epidermal thickness, as well as a 
thinning of the granular epidermal layer compared with WT 
mice, more pronounced in the foot pads versus dorsal skin 
(Fig. 5, A and B). In addition, Satbl - ' - mice showed signifi- 
cantly (P < 0.05) decreased epidermal proliferation compared 
with WT controls (Fig. 5 C). 

To gain mechanistic insight into the role of Satbl in the 
regulation of the keratinocyte terminal differentiation, samples 
of total RNA isolated from the skin of newborn Satbl - ' - and 
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Figure 3. Satbl is a direct target of p63 in keratinocytes. Primary kera- 
tinocytes isolated from El 6.5 embryos or newborn mice were processed 
for ChIP analysis with an antibody against the p63 protein or purified 
rabbit IgGs. HaCaT cells were used in cotransfection experiments with 
Satbl promoter-driven reporter construct. (A and B) Quantitative RT-PCR 
analysis of two distinct regions of the Satbl (a predicted high-affinity p63- 
binding site and a negative control site) showing a specific high-affinity 
p63-binding site at the Satbl promoter. In A, the position 1 refers to 5' of 
Satbl transcript (AK037740). The uppercase letters of the sequence show 
the p63 putative binding sites in the promoter region of the Satbl gene 
chosen for the quantitative PCR analysis after ChIP. The promoter region 
of the Cldn 1 was used as a positive control. The input levels of unprecipi- 
tated chromatin DNA were used as loading controls. Error bars represent 
SEM, and three independent experiments were run in triplicate. (C] HaCaT 
keratinocytes were cotransfected with the luciferase reporter plasmid con- 
taining a mouse Satbl promoter fragment and pANp63 expression plas- 
mid or control pcDNA3 plasmid. The increase in the pSatbl-luc activities 
by cotransfection with pANp63 compared with the control pcDNA3 was 
~3.3 fold (mean ± SEM, n = 3; *, P < 0.05, Student's ftest). 
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Figure 4. Alterations in the conformation of the 5 Mb chromatin domain containing EDC in the epidermis of p63 -/ ~ and Satbl mice. The skin of E 1 6.5 
p63 -/, ~, Satbl and corresponding WT mice was processed for 3D FISH analyses, which were correlated with changes in gene expression determined 
by quantitative RT-PCR. (A) A schematic structure of the 5 Mb domain on mouse chromosome 3 (mChr3) containing the EDC locus, Rps27, and Gabpb2 
genes. 3D FISH DNA probes detecting the corresponding domains are shown in green [Loricrin], pink [Rps27], and yellow (Gabpb2). Genes selected for 
quantitative RT-PCR analysis (shown in E) are shown in red. LEP, late-cornified envelope protein. (B) Multicolor 3D FISH with BACs covering the Rps27, Lor, 
and Gabpb2 in the epidermal cells of p63~/~, Satbl and corresponding WT mice at El 6.5 (representative single Z sections). Bar, 2 urn. (C) 3D FISH 
distances between the Rps27and Lor normalized to the radius of each nuclei in basal epidermal cells of p63 -/ ~, Satbl and corresponding WT mice. 
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Figure 5. Satbl knockout mice show alterations in the epidermal structure and Loricrin expression. Cryosections of the footpad skin of newborn Satbl 
and WT mice were processed for morphometric and immunofluorescent analyses. (A) Alterations in the structure and thinning of granular layer (arrows) in 
the epidermis of Satbl mice compared with WT mice. Bars, 25 urn. (B) Significant (*, P < 0.05) decrease of the epidermal thickness in the epidermis 
of Satbl -/, ~ mice compared with WT mice. Measurements of the epidermal thickness were performed in three Satbl _ '~ and three WT mice. 20 measure- 
ments were performed in each mouse. (C) Significant (*, P < 0.05) decrease of cell proliferation in the epidermis of Satbl ~'~ mice versus WT controls. 
The percentages of the Ki-67+ nuclei were determined in the epidermis of two Satbl ~'~ and two WT mice; 40-50 nuclei were analyzed in each animal. 
(B and C) Error bars represent SEM. (D) Decrease of Loricrin expression in the epidermis of Satbl versus WT mice. Bars, 25 urn. (E) Agilent microarray 
data demonstrating changes in gene expression in the EDC between P0.5 Satbl and WT mice. Asterisks indicate the fold changes in gene expression 
levels validated by quantitative RT-PCR. 



WT mice were processed for microarray analysis (Sharov et al., 
2006). In addition, the formaldehyde-cross-linked chromatin 
fragments from primary epidermal keratinocytes were immuno- 
precipitated with anti-Satbl antibody, and a NimbleGen MM8 
Deluxe Promoter HX1 array analysis was performed. ChlP- 
on-chip analysis showed that Satbl does not always bind in the 
vicinity of the transcription start sites of the genes. In some 
cases, it binds sites remote from the gene proximal regulatory re- 
gions (Tables S4-S6), which is consistent with its role as a ge- 
nome organizer important for the proper higher-order chromatin 



folding previously shown (Cai et al., 2003, 2006; Wang et al., 
2009). Microarray and ChlP-on-chip data were merged, and 
2,503 genes whose expression was significantly changed in 
Satbl~'~ mice (>1. 8-fold) and showed Satbl binding within 
200 kb from the corresponding transcription start sites were 
selected as targets for Satbl in keratinocytes (Fig. S3 A and 
Tables S4 and S5). 

Interestingly, bioinformatic analyses on the Satbl target 
genes involved in many biological processes (viz., regulation 
in cell adhesion/ECM remodeling, metabolism, cytoskeleton, 



Pairwise comparisons represent a significant increase in the Lor-Rps27 distances between the E 1 6.5 WT versus p63~^~ or Satbl ~'~ mice (mean ± SEM, 
n = 60; *, P < 0.01 , Newman-Keuls test after a one-way ANOVA test). (D) 3D FISH distances between the Rps27and Lor and the Lor and Gabpb2 in basal 
epidermal cells and dermal cells of p63~/~, Satbl and corresponding WT mice. Pairwise comparisons represent a significant increase in the Lor-Rps27 
distances between the E 1 6.5 WT versus p63~ /_ or Satbl mice (mean ± SEM, n = 60; P < 0.01 , Newman-Keuls test after a one-way ANOVA test). 
(E) Quantitative RT-PCR analyses of gene expression in the 5 Mb domain containing EDC in the epidermis of E 1 6.5 p63~^~ and Satbl mice normalized 
to the expression levels in the corresponding WT mice. White slashes represent the interruption of these bars to indicate that they show the fold changes in 
gene expression levels higher than the maximum value in the y axis. The exact fold changes are indicated under each bar. SEM is shown in the appropriate 
scale (relative to the y axis). Three independent experiments were run in triplicate. 
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Figure 6. Treatment with Satbl -expressing lentivirus partially rescues alterations of epidermal phenotype in skin explants caused by p63 deficiency. The 

skin of p63 +/ ~ and WT El 3.5 embryos was processed for modulation of gene expression in organ culture with p63 siRNA and Satbl -expressing lentivirus. 
(A and B) Morphology of E 1 3.5 skin explants of p63 +/ ~ mice treated with combinations of p63 or control siRNAs and Satbl -expressing or control lentiviruses 
(LV) for 48 h. The decrease of the epidermal thickness induced by p63 siRNA is significantly (**, P < 0.01 ) rescued by the cotreatment with Satbl -expressing 
lentivirus compared with corresponding controls. Epidermal thickness was measured in three samples for each experimental group; 60 measurements were 
conducted in each sample. (C and D] Cotreatment with Satbl -expressing lentivirus prevents the decrease of cell proliferation in the epidermis of skin explants 
induced by p63 siRNA compared with controls (*, P < 0.05). The percentages of Ki-67+ cells were analyzed in three samples for each experimental group; 
30 cells were analyzed in each sample. (A and D) The dotted lines show the position of the basement membrane separating the epidermis and dermis. 
(E and F) Quantitative RT-PCR for Loricrin and Satbl (E) and immunostaining for Loricrin (F). Levels of Loricrin transcripts and the expression of protein de- 
creased after p63 siRNA treatment restored in the epidermis after cotreatment with Satbl -expressing lentivirus compared with the controls (F, arrows). Treat- 
ment with Satbl -expressing lentivirus resulted in a marked increase of Satbl transcript levels in the epidermis after decrease induced by p63 siRNA (E). 
(A, D, and F) The epidermal-dermal junction is shown by dotted lines. (B, C, and E) Error bars represent SEM. Bars: (A and F) 50 pm; (D) 25 \im. 



signaling, and transcription) revealed that a large number of those 
genes are grouped into large genomic loci (>0.5 Mbp each) of 
gene clusters, containing >20 functionally relevant genes per 
locus, and that such genomic loci include epidermis-specific loci 
(Fig. S3 B and Tables S4 and S5). In fact, analyses of 33 such loci 
that are distributed genome wide showed significant preferential 
binding (P < 0.05) of Satbl and Satbl -dependent changes in gene 
expression in the keratinocyte-specific loci (EDC, keratin type 
I/II loci, and keratin-associated protein locus) located on mouse 
chromosomes 3, 11, 15, and 16 (Fig. S3 B). Many of those genes 
found in the keratinocyte-specific loci were altered in expression 
upon p63 ablation (Fig. 5 E and Tables SI and S2). 

Microarray analysis revealed a profound decrease in gene 
expression in all keratinocyte-specific loci in Satbl - ' - mice 
compared with WT controls (Fig. 3 B and Tables S4 and S5). 
Similarly to p63 - ' - mice, the expression of KrtlO, Krt23, and 
Krt78 decreased in Satbl - ' - mice compared with WT controls 
(Tables SI and S4). In the EDC locus, Satbl ablation resulted in 
decreased expression of a large number of genes that belong to 



the small proline-rich and late-cornified envelope families 
involved in the control of terminal keratinocyte differentiation, 
whereas an increase in expression was detected for selected 
members of the S100 family (Fig. 5 E). Consistent with micro- 
array and quantitative RT-PCR data (Fig. 4 E), immunohisto- 
chemical analysis revealed markedly reduced expression of the 
Loricrin in Satbl-null epidermis versus the controls (Fig. 5 D). 
Thus, alterations of epidermal morphology in Satbl - ' - mice 
and marked similarity in the Satbl - and p63-dependent changes 
in gene expression in four major keratinocyte-specific loci 
suggest that Satbl is indeed a significant component of p63- 
regulated molecular network in the epidermis. 

Gain of Satbl expression partially rescues 
the epidermal phenotype in embryonic skin 
explants of p63-def icient mice 

To assess whether the restoration of Satbl expression is capable 
of partially rescuing the skin phenotype of p63 - ' - mice, El 3. 5 
skin explants isolated from heterozygous (+/— ) p63 knockout 
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mice were cultured ex vivo and treated with p63 siRNA in com- 
bination with either Satbl -expressing or control lentiviruses 
(Fig. 6). Treatment of p63 +7- skin explants with p63 siRNA re- 
sulted in an ^80-90% decrease of the p63 mRNA levels com- 
pared with WT skin (Fig. S2 A). p63 siRNA treatment also 
caused marked decrease of the Satbl mRNA levels, which, 
however, were restored by cotreatment with Satbl -expressing 
but not the control lentiviruses (Fig. 6 E). 

In contrast to the p63 +/ ~ skin explants treated by control 
siRNA, treatment with p63 siRNA alone or in combination with 
control lentiviras caused a marked decrease of the epidermal 
thickness, cell proliferation, and Loricrin expression, thus repro- 
ducing in part the skin phenotype of p63 -7- mice (Fig. 6, A-F). 
However, cotreatment of p63 +/ ~ skin samples with p63 siRNA 
and Satbl-expressing lentiviras resulted in a significant increase 
in epidermal thickness, cell proliferation, and Loricrin expression 
compared with the samples treated with p63 siRNA and control 
lentiviras (Fig. 6 A-F). Restoration of epidermal proliferation 
after treatment with Satbl-expressing lentiviras was consistent 
with data showing the reduced cell proliferation in the epidermis 
of Satbl - ' - mice compared with WT controls (Fig. 5 C). There- 
fore, it is likely that Satbl plays a role in regulating both kera- 
tinocyte proliferation and differentiation and is capable of partially 
restoring the epidermal phenotype of p63-deficient mice. 

Satbl binds to the distinct genomic 
regions within the EDC locus 
in keratinocytes 

Because the distance between the Lor and Rps27 in the 5 Mbp 
chromatin domain is increased in Satbl -7- epidermis in com- 
parison with WT mice (Fig. 4, B-D), we hypothesized that 
Satbl may bind to specific sites in the EDC locus to regulate the 
EDC chromatin conformation. To determine whether any re- 
gions within the EDC locus are direct binding targets of Satbl, 
the formaldehyde-cross-linked chromatin fragments from kera- 
tinocytes isolated from E16.5 embryos or newborn skin were 
immunoprecipitated with anti-Satbl antibody, and ChlP-on-chip 
analysis was performed (Fig. 7 A). 

Analyses of the ChlP-on-chip data revealed Satbl binding 
to many genomic sites, and these sites are found enriched in 
three distinct regions within the EDC locus (at the 5' end do- 
main adjacent to the gene desert region, the central domain, and 
the 3' end domain). Importantly, Satbl-binding sites in the cen- 
tral domain correspond to the small proline-rich protein (Sprr) 
and late-cornified envelope (Lee) gene families (Fig. 7 A and 
Table S6), which are activated during terminal keratinocyte dif- 
ferentiation (Marshall et al., 2001; Martin et al., 2004; Bazzi 
et al., 2007; Brown et al., 2007). ChlP-on-chip data were further 
validated by ChlP-quantitative PCR, and the results confirmed 
that, in both E16.5 and newborn mice, Satbl binds to the proxi- 
mal regulatory regions of the Sprrlb and Lce3b genes located at 
the central EDC domain (Fig. S4). 

Satbl compresses chromatin conformation 
of the EDC locus and its central domain 

To further examine the potential roles of Satbl in chromatin re- 
modeling and the execution of gene expression programs in 



epidermal progenitor cells, 3D FISH analysis of the chromatin 
conformation was performed in the epidermis of E16.5 and 
P0.5 Satbl -7- and age-matched WT mice, focusing on total 
EDC and its central domain containing a large number of the 
keratinocyte-specific genes (such as genes belonging to the 
Sprr and Lee families as well as Lor and Inv). Using DNA 
probes corresponding to S100a6 (5' end), Lor-Lce3c (the 
central domain), and SlOOalO (3' end), we analyzed the 3D 
structure of the EDC locus by measuring the volume of its 
central domain as well as distances between its 5 ' end, central 
domain, and 3' end. We found that in E16.5 or P0.5 Satbl -7- 
and WT mice, the EDC locus was localized in the nuclear 
interior of epidermal progenitor cells, and Satbl ablation did 
not significantly alter its radial positions compared with WT 
mice (Fig. S5, A and B). 

However, the most striking feature of the EDC confor- 
mation in the El 6.5 and P0.5 Satbl -7- epidermis, in com- 
parison with age-matched WT mice, was the significantly 
(P < 0.001) expanded central domain, determined as the vol- 
ume of the Lor-Lce3c signal (representative Z optical sections 
of the corresponding DNA FISH images are shown in Fig. 7 B, 
and continuous Z section images are shown in Fig. S5 E). As 
evident from the distribution of the data, we consistently ob- 
served the increased volume of the central domain in each in- 
dividual nuclei of epidermal cells in P0.5 Satbl -7- mice 
compared with aged-matched WT mice (Fig. 7 C). In contrast 
to the epidermal cells, a lack of changes in the volume of the 
central EDC domain was seen between dermal cells of 
Satbl -7- and WT mice (Fig. S5 D), thus suggesting that the 
effects of the Satbl ablation on its conformation are indeed 
keratinocyte specific. 

In addition, we found an overall expansion in the length of 
the whole EDC in P0.5 Satbl -7- epidermis (P < 0.01) compared 
with age-matched WT mice (Figs. 7 [C and D] and S5 C), which 
is likely a result of elongation of its 5' domain, determined as a 
distance between the S100a6 and the center of the Lor-Lce3c 
signal (Table S7). These results are consistent with data demon- 
strating the significant elongation of the distance between Lor 
and Rps27 genes in Satbl - ' - embryos versus WT mice (Fig. 4, 
B-D) and suggest that Satbl ablation leads to significant de- 
compression of the locus. 

Interestingly, Satbl binding at specific sites within the 
EDC locus was not always detected at the proximal promoter 
regions of the genes whose expression was markedly altered 
in the Satb -7- epidermis (e.g., Lor, Fig. 7 A and Table S6). 
Our 3D FISH data indicate that Satbl compresses the chro- 
matin conformation at the EDC locus, including the regions 
enriched in genes activated during terminal keratinocyte dif- 
ferentiation (Fig. 7, B-D). We observed marked differences 
in the expression of the genes located in the EDC central do- 
main upon either Satbl or p63 ablation (Figs. 4 E and 5 E). 
Therefore, it is likely that Satbl regulates gene expression in 
the EDC locus in epidermal progenitor cells by establishing 
its tissue-specific chromatin conformations and potentially by 
forming a densely looped chromatin structure, rather than by 
directly targeting proximal cis-regulatory elements of indi- 
vidual genes. 



pB3. Satbl , and tissue-specific chromatin remodeling • Fessing et al. B33 



o 
_l 

"D 

0> 

o 



chr3 EDC 


1 1 


5-end Central domain 3 -end 
i 1 i 1 i — i 


f 


i i 

incl. S100a6) 


■ iii 

(Lor-Lce3c) (incl. S100a10) 
















































Gene "desert" 
























i 








—1 L 








1 1 1 










1 





B 



Rps27 S100 
Family 



Lor SPRR 
Family 



LEP 

Family 



FN S100 
Family 



Gabpb2 



WT 



Satbl-/- 

































S100a1l 


I 





Total EDC length (P0.5) 



Central domain (P0.5) 

Satbl-/- 
WT 




20 30 40 

Loci 



WT 



r~ 

5' end 



S100a6 



Satbl -/- 



S100a6 



E16.5 

1113 ±84.6 nm 



center 



WT 



3 end 



5 end 



Lor Lce3c S100a10 

1211 ± 112 nm 



S100a6 




Satbl -/- 

r - 



Lor Lce3c S100a10 

**p<0.001 



S100a6 



P0.5 

1046 ± 88.4 nm 



I 



center 



1 

3 end 



-^53±4.6 x10^n^ ^ 0 - 



Lor 



Lce3c S100a10 



1281 ± 104.9 nm' 
1 




Lor Lce3c S 100a 10 

* p<0.01 ** p<0.001 



Figure 7. Satbl binds the central EDC domain and regulates its conformation in epidermal cells. Primary mouse keratinocytes were processed for ChlP- 
on-chip analyses with Satbl antibody, and skin cryosections of E 1 6.5 and P0.5 Satbl and WT mice were processed for 3D FISH analyses of the EDC 
chromatin structure. (A) ChlP-on-chip analysis of the Satbl binding to the distinct genomic regions of the 5 Mb chromatin domain on mouse chromosome 
3 (chr3) containing EDC in primary keratinocytes. Enrichment of the Satbl -binding sites at the distinct EDC regions including the central domain. Lack of 
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Discussion 

Here, we show that p63 transcription factor governs a program 
of the development of stratified epithelia (Koster and Roop, 
2007; Truong and Khavari, 2007; Crum and McKeon, 2010) via 
regulating expression of the genes that control higher-order 
chromatin remodeling, covalent histone modifications, and nu- 
clear assembly. Among the genes that regulate chromatin struc- 
ture and gene expression, we focused on Satbl, which is known 
to regulate large-scale chromatin remodeling in selected cell 
types and progenitor cells (Cai et al., 2003, 2006; Kumar et al., 
2007; Han et al., 2008; Agrelo et al., 2009). Satbl was identi- 
fied among the genes involved in chromatin remodeling that 
were down-regulated in p63~'~ epidermis, consistent with data 
previously reported using cultured epithelial cells (Carroll et al., 
2006; Truong et al., 2006; Pozzi et al., 2009). Our data unveil a 
novel mechanism by which p63 functions in epidermal develop- 
ment, at least in part, via direct regulation of Satbl . 

We demonstrate that in keratinocytes, Satbl shows prefer- 
ential binding to four major tissue-specific gene loci (EDC, kera- 
tin type l/ll loci, and keratin-associated protein locus). Satbl 
ablation, in striking similarity to p63 deficiency, results in pro- 
found alterations of expression of a large number of epidermis- 
specific genes. Furthermore, we show that the restoration of Satbl 
levels partially rescues the epidermal phenotype of embryonic 
skin explants from p63-deficient mice, further suggesting Satbl 
as an important part of the p63-regulated network in the skin. 

Among many chromatin-remodeling factors (such as 
Smarca4/Brgl, Mi-2a, Mi-2(3, etc.) that are altered in p63" /_ 
mice and other p63 targets previously identified (Koster and 
Roop, 2007; Vigano and Mantovani, 2007), Satbl represents an 
important functional component of the p63-regulated genetic 
network because of its ability to reorganize chromatin structure 
and regulate higher-order chromatin remodeling in epidermal 
progenitor cells. In fact, our expression data together with ChlP- 
on-chip data for Satbl -binding sites indicate that once Satbl is 
induced by p63, it establishes its own regulatory network, bind- 
ing to the target gene loci and controlling the expression of 
^2,500 genes, including those not associated with the epider- 
mis-specific gene loci. Among these genes, cell cycle-associated 
genes are of interest because we found that Satbl has an effect 
in promoting cell proliferation in the developing epidermis, 
which is consistent with data obtained from the studies on ag- 
gressive breast cancer cells (Han et al., 2008). The mechanisms 
underlying Satbl -induced cell proliferation need to be studied 
in the future. 

Differences in gene expression profiles in keratinocytes 
observed after either p63 or Satbl ablation in many gene groups 



may suggest that the roles of p63 and Satbl in the control of ex- 
pression of the genes located outside of the epidermis-specific 
loci are quite distinct. However, we must be aware that p63 abla- 
tion, which deregulates many key transcriptional repressors/ 
activators and chromatin-remodeling factors, in addition to Satb 1 , 
could alter expression of many Satbl target genes beyond the 
effect of Satbl ablation alone. This could explain, at least for some 
gene groups, why their expression profiles are quite distinct be- 
tween p63~'~ and Satbl mice. Nevertheless, we provide com- 
pelling evidence that within a spiderweb of a gene regulatory 
network governed by p63, Satbl, as its direct downstream target 
gene, constitutes an important component driving epidermal de- 
velopment via regulating the epidermis-specific gene loci. It still 
remains to be investigated whether Satbl acts together with p63 
or whether it acts independently from p63 in controlling the 
establishment of specific conformation of the EDC and other 
genomic loci in keratinocytes. 

We show here that Satbl binds to distinct genomic regions 
including the central domain of the EDC locus, which consists of 
a large number of genes activated during terminal keratinocyte 
differentiation (Martin et al., 2004; Brown et al., 2007) and which 
regulates the establishment of its specific conformations in kera- 
tinocytes. It has been previously shown that intranuclear EDC 
positioning in cultured keratinocytes is cell type specific, whereas 
their treatment with 5-azacytidine and sodium butyrate results in 
an increase of gene expression within the EDC (Elder and Zhao, 
2002; Williams et al., 2002). These and other studies suggest that 
changes in local epigenomic modification of the EDC gene locus 
are required for efficient expression of keratinocyte-specific 
genes (Frye et al., 2007; Ezhkova et al., 2009; Sen et al., 2010). 

Furthermore, we demonstrate that Satbl deficiency results 
in decompression of the EDC locus, specifically its central do- 
main. These alterations are associated with a marked decrease 
in expression of terminal differentiation-associated genes and 
are accompanied by morphological changes in the epidermis 
(thinning of the granular layer and decrease of the epidermal 
thickness). Lack of similar dynamics in the nuclei of dermal 
cells suggests that the developmentally regulated remodeling of 
the higher-order chromatin structure in the EDC associated with 
increase of its transcription activity is indeed a tissue-specific 
event and that other factors may control EDC condensation in 
dermal cells. 

During the execution of cell differentiation programs, tissue- 
specific genes organized into large genomic loci are expressed in a 
tightly coordinated manner to achieve the proper balance of the 
distinct proteins in differentiated cells (Kosak and Groudine, 
2004). Increased evidence of data suggests that for such coordi- 
nated regulation, distinct spatial organization of tissue-specific 



Satbl binding to the gene desert region between the 5' end and the central domain of the locus. LEP, late-cornified envelope protein. (B] 3D FISH analyses 
of the volume of central EDC domain as well as the EDC length with probes detecting the 5' and 3' ends of the EDC [S100a6 and S 100a 10 genes, respec- 
tively) and the central domain (Lor-Lce3cj in the basal epidermal cells of Satbl and WT mice. Representative single Z optical sections of newborn WT 
and Satbl knockout keratinocyte nuclei are shown. Consecutive Z sections of the same nuclei are shown in Fig. S5 E. Arrows indicate the FISH signals. Bar, 
2 urn. (C) The volume of the EDC central domain and the length of the total EDC locus were measured and compared between nuclei of basal epidermal 
cells of Satbl and WT mice. The distribution of the data after measuring the total EDC length with or without normalization to the corresponding nuclear 
radii (left and middle, respectively) and the volume of the central domain (right) in 60-70 loci of basal epidermal cells at P0.5 Satbl and WT mice is 
shown. Individual loci are shown in the x axis. (D) Statistical analysis shows a significant increase of the volume of central EDC domain in Satbl mice 
compared with that of WT mice. The length of the total EDC increases in P0.5 Satbl mice compared with WT mice. 
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gene loci is required, and Satbl was shown to be an essential fac- 
tor in controlling formation of the specific loop structures in many 
tissue-specific gene loci (Cai et al., 2003, 2006). It has been 
shown that during T H 2 lymphocyte activation, Satbl promotes 
the formation of a specific 3D structure of the cytokine locus, in 
which chromatin is folded into numerous small loops that bring 
proximal and distal gene regulatory elements together to activate 
the gene expression (Cai et al., 2006). More recent data demon- 
strate that Satbl also controls the formation of the chromatin 
loops in the (3-globin locus (Wang et al., 2009). We hypothesize 
that Satbl may form a similar chromatin structure, characterized 
by dense chromatin loops, in the EDC central domain and the 
flanking regions, which may result in reduction of the total physi- 
cal volume of chromatin in these regions. Such a higher-order 
transcriptionally active chromatin structure may facilitate the 
accessibility of protein complexes to target gene loci as well as 
interactions between proximal and distal regulatory gene regions. 

If this hypothesis is correct, in the absence of Satbl, the 
EDC locus cannot be properly folded, and, therefore, this region 
on mouse chromosome 3 will remain more decompressed (or 
elongated) compared with WT cells and will not be able to sup- 
port efficient expression of the terminal differentiation-associated 
genes. To test this hypothesis, additional studies using the chro- 
matin conformation capture techniques will be required to iden- 
tify the role of Satbl in the formation of such chromatin loop-like 
structures within the EDC. However, significant changes in the 
volume of the central EDC domain and total EDC length upon 
Satbl ablation provide strong support for this hypothesis. 

Importantly, these data are also consistent with previous 
observations that strong Satbl binding may not be necessarily 
detected at the proximal promoter regions of its target genes. This 
is because the main binding targets of Satbl are the specialized 
DNA domains called base unpairing regions (Kohwi-Shigematsu 
and Kohwi, 1990; Bode et al, 1992), which are enriched in the 
gene-rich regions throughout the genome (Cai et al., 2006). For 
instance, Satbl binding was undetected at the proximal regula- 
tory region of the Lor gene per se, and yet Lor appears to be regu- 
lated in the context of Satbl -dependent 3D chromatin folding. 

The data presented in this paper link, for the first time, 
the functions of a master regulator of tissue development (viz., 
p63) with a protein that controls higher-order chromatin re- 
modeling (viz., Satbl), and this link provides an important 
foundation for further studies of the role of p63 as a potential 
upstream regulator of Satbl activity in other cell types. Be- 
cause the Satbl-dependent step in the p63-mediated program 
of higher-order chromatin remodeling in the EDC locus is 
highly important for the proper execution of differentiation 
programs in epidermal progenitor cells, it remains to be deter- 
mined whether similar mechanisms underlie the involvement 
of p63 and Satbl during the developmentally regulated activa- 
tion of gene expression in other keratinocyte-specific loci as 
well as during activation of the epidermis-specific gene ex- 
pression programs in epithelial stem cells in adult skin (Ito 
et al., 2007; Watt et al., 2008; Blanpain and Fuchs, 2009). 
Importantly, recent data indeed revealed that Satbl has a criti- 
cal role in regulating the differentiation potential of embry- 
onic stem cells (Savarese et al., 2009). 



Collectively, our data demonstrate that Satbl is an important 
hitherto unrecognized component of the global regulatory network 
of p63, a master regulator of epidermal morphogenesis. These 
data provide a novel fundamental mechanism showing how mas- 
ter regulators of tissue morphogenesis establish tissue- and devel- 
opmental stage-specific patterns of chromatin organization and 
remodeling to control the fate of multipotent progenitor cells. 

Materials and methods 

Experimental animals and tissue collection 

All animal studies were performed under protocols approved by the Boston 
University and University of California Berkeley Institutional Animal Care 
and Use Committee and the Home Office Project License. C57BI/6 mice 
were purchased from Charles River. Skin samples were collected from mice 
at distinct days of embryonic and postnatal development (E 1 6.5 and P0.5). 
p63~/~ and p63 +/+ control embryos at El 6.5 were obtained by breeding 
p63 +/ ~ animals purchased from The Jackson Laboratory. Genotyping 
of mice was performed using PCR, as recommended by the supplier. 
For each developmental stage, six to seven samples were collected. Tis- 
sue samples were frozen in liquid nitrogen and processed for immuno- 
fluorescent or microarray analyses as previously described (Sharov et al., 
2003, 2005b). 

In brief, for FISH and immunofluorescent analysis of 3D preserved 
nuclei, skin samples were processed through the fixation protocol with form- 
aldehyde and gradual equilibration in buffered sucrose solutions according 
to previously published recommendations (Solovei etal., 2009). The embryos 
or tissues were fixed in 4% formaldehyde at 4°C overnight and incubated in 
50 mM ammonium chloride solution at room temperature for 5 min followed 
by washes with 0.1 M phosphate buffer (pH 7.0) and by incubations in 
5 and 1 2.5% sucrose solution in the phosphate buffer (1 h at room tempera- 
ture) and in 20% sucrose solution (overnight at 4°C). Finally, the samples 
were frozen and embedded into optimal cutting temperature medium. 

LCM, microarray, and quantitative RT-PCR analysis 

LCM of whole mouse epidermis at selected stages of development was per- 
formed followed by RNA isolation, amplification, and microarray analysis 
as previously described (Sharov et al., 2006). In brief, 8-pm-thick frozen 
skin cryosections were dehydrated to preserve RNA integrity, and LCM was 
performed from the epidermis of WT and p63~^~ mice at El 6.5 using an 
LCM microscope (Life Technologies). Total RNAs were isolated using the 
PicoPure RNA isolation kit (Life Technolgoies) followed by two rounds of lin- 
ear RNA amplification using the RiboAmp RNA amplification kit (Life Tech- 
nologies). Equal amounts of RNA from each sample were labeled by Cy3 
using a fluorescent labeling kit (Agilent Technologies). Qualify and size dis- 
tribution of targets was determined by RNA 6000 Nano Laboratory- 
on-chip assay (Agilent Technologies), and quantification was determined 
using a microscale spectrophotometer (NanoDrop). After one round of lin- 
ear amplification in all analyses, Universal Mouse Reference RNA (Agilent 
Technologies) was used as a control. All microarray analyses were per- 
formed by MOgene, LC using the 4 1 K Whole Mouse Genome 60-mer Oligo 
Microarray kit (Agilent Technologies). Real-time PCR of unamplified refer- 
ence RNA and reference RNA obtained after two rounds of amplification 
was used for validation of possible alterations in gene expression caused by 
the amplification procedure. All microarray data on gene expression were 
normalized to the corresponding data obtained from the reference RNA. 
Two independent datasets were obtained from WT and transgenic mice, and 
p values were calculated by Feature Extraction software (version 7.5; Agilent 
Technologies) using distribution of the background intensity values to signal 
intensity and using a Student's f test. For quantitative RT-PCR analysis of RNA 
isolated and amplified after LCM, PCR primers were designed on Beacon 
Designer software (Table S6; PREMIER Biosoft International). Real-time PCR 
was performed using iQ SYBR green Supermix and the MyiQ Single-Color 
Real-Time PCR Detection System (Bio-Rad Laboratories). Differences between 
samples and controls were calculated using the Gene Expression Macro pro- 
gram (Bio-Rad Laboratories) based on the AAC, equitation method with glyc- 
eraldehyde 3-phosphate dehydrogenase and peptidylprolyl isomerase A 
genes as internal housekeeping controls. 

3D FISH and immunofluorescent analyses 

Labeled DNA probes were prepared, and 3D FISH analysis was performed 
on 20-um tissue sections as previously described (Solovei et al., 2009). 
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Individual or pooled bacterial artificial chromosome (BAC) DNA for the 
genetic loci of interest was used for the probe synthesis (Table S3). The 
probes were synthesized by nick translation using Biotin-dUTP, FITC-dUTP, 
or Cy3-dUTP. Whole chromosome paints were labeled with Biotin-dUTP by 
degenerate oligonucleotide-primed PCR. The frozen sections were briefly 
dried, heated in 10 mM sodium citrate (pH 6.0) to partially reverse cross- 
linking, and equilibrated in 50% formamide/2x SSC. Labeled DNA probe 
mixtures were dissolved in hybridization buffer and applied on the tissue 
under glass chambers. Genomic and probe DNA were denatured by incu- 
bating slides at 85°C for 5 min, and hybridization was performed at 37°C 
for 2 d. After hybridization, slides were washed three times at 2x SSC for 
1 0 min at 37°C followed by one wash in 0. 1 x SSC for 1 0 min at 60°C. 
When biotinylated probes were used, the slides were incubated with avidin- 
Alexa Fluor 488 (Invitrogen), streptavidin-Cy3, or streptavidin-Cy5 (Rockland 
Immunochemicals). DNA was stained with DAPI (Sigma-Aldrich), and 
slides were embedded using VectaShield medium (Vector Laboratories). 
Immunofluorescent analysis was performed on 1 0-20-pm quickly frozen or 
3D preserved tissue with antibodies (Table S10) using standard protocols 
(Fessing etal., 2006). 

Microscopy and image analyses 

For analysis of epidermal morphology, images of the tissue after staining 
with hemotoxylin were acquired at room temperature using a microscope 
(Eclipse 50i; Nikon) equipped with a Plan Fluor 20x/0.50 or 40x/0.75 
objective lens (Nikon), VisiCam 3.0 camera (VWR International), and Visi- 
Cam Image Analyzer software. The morphometric analysis was performed 
using ImageJ software (National Institutes of Health). Immunofluorescent 
images were acquired at room temperature using a microscope (Eclipse 
50i), EXi Aqua camera (Qlmaging), and Image-Pro Express software (ver- 
sion 6.3; Media Cybernetics). 

3D FISH images were collected using a laser-scanning confocal 
microscope (510 META; Carl Zeiss), Plan Apochromat 63x/l .4 oil differ- 
ential interference contrast objective lens, and laser-scanning microscopy 
software (510 META). Imaging was performed at room temperature using 
Immersol immersion oil (Carl Zeiss). For 3D image analysis, the nuclei were 
scanned with a Z axial distance of 200 nm yielding separate stacks of 
8-bit grayscale images for each fluorescent channel with a pixel size of 
1 00-200 nm. For each optical section, images were collected sequentially 
for all fluorophores, and axial chromatic shift was corrected for each chan- 
nel using 500-nm TetraSpeck beads (Invitrogen; Solovei et al., 2009). The 
images were processed and analyzed using ImageJ software. For the illus- 
tration purposes, 3D FISH and immunofluorescent images of the skin cryo- 
sections were adjusted using the levels and brightness/contrast tools in 
Photoshop (CS4; Adobe); separate adjustments were applied to every 
pixel in each RGB channel. 

To determine radial positions of the loci, coordinates of at least 200 
random points from the nuclear surface depicted by DAPI counterstain were 
acquired and exported to Excel (Microsoft), where the nuclear geometric 
center was calculated by coordinate averaging (Hofers et al., 1993). The 
mean nuclear radius for each nucleus was calculated as a mean of dis- 
tances from the nuclear geometric center to all surface points. Distances be- 
tween the centers of the nucleus and distinct loci were calculated on the 
basis of their spatial coordinates and were normalized to the mean nuclear 
radius using Excel (Ronneberger etal., 2008). Distances between individual 
gene loci were calculated after determination of the spatial coordinates of 
the corresponding centers of the signals using ImageJ software. The dis- 
tances were normalized to the mean radii of the corresponding nuclei, de- 
termined as described above. 

The volume of the central EDC domain was calculated using the ap- 
proach previously described (Rauch et al., 2008). First, the threshold of the 
fluorescent signal was determined via selection of 1 % of the brightest vox- 
els from 3D intensity histograms using ImageJ software. Second, using this 
threshold, the number of voxels within each fluorescent signal was calcu- 
lated. Next, the physical volume indicated by the fluorescent signal was 
determined as a sum of the volumes from all voxels after correction for opti- 
cal aberrations by normalizing the data to the volume of beads of distinct 
size (100 nm, 500 nm, and 4 urn; TetraSpeck). 

Embryonic tissue culture 

Whole dorsal skin culture from El 3.5 p63 +/ ~ embryos were prepared and 
maintained as previously described (Botchkarev et al., 1 999). In brief, skin 
samples were dissected from embryos and cultured in the 6-well plates 
containing Williams' E culture medium for 48 h. Within 2 h after the culture 
preparation, tissue samples were transfected with Stealth siRNA, recogniz- 
ing all isoforms of p63 (Antonini et al., 2006), or medium containing 



negative control Stealth siRNA (Invitrogen) at 100 nM final concentration 
using the RNAiMax reagent (Invitrogen). Samples were also incubated 
with combinations of p63 siRNA or control siRNA and Satbl -expressing 
lentivirus or control lentivirus (10 ug/ml) for 48 h. The tissue was frozen 
and embedded into optimal cutting temperature medium as described 
in the Experimental animals and tissue collection section for subsequent 
immunofluorescent analyses. The penetration of the RNAi into the tissue 
under the used experimental conditions was demonstrated using BLOCK-iT 
Fluorescent Oligo (Invitrogen) in parallel experiments. 

Western blot analysis 

Proteins were extracted from snap-frozen skin samples or cultured keratino- 
cytes with lysis buffer as previously described (Antonini et al., 2006). 5 pg 
of protein was processed for Western blot analysis followed by incubation 
with primary antibodies against p63 or Erk (Table S10) overnight at 4°C. 
HRP-tagged IgG antibody was used as a secondary antibody (1:5,000; 
Thermo Fisher Scientific). Antibody binding was visualized using an ECL 
system (SuperSignal West Pico kit; Thermo Fisher Scientific) followed by 
autoradiography with x-ray film (Thermo Fisher Scientific). 

ChIP and ChlP-on-chip assays 

ChIP assay was performed using epidermal keratinocytes isolated 
from El 6.5 or newborn mouse skin with anti-p63 or anti-Satbl antibodies 
(Table S8) as previously described (Antonini et al., 2006; Nguyen et al., 
2006). PCR analysis for enrichment in the promoter sequences of Satbl, 
Cldnl, Lor, Sprrlb, Lce3b, or Inv genes in precipitated DNA was per- 
formed using the PCR primers listed in Table S7. For ChlP-on-chip analysis, 
precipitated DNA was analyzed by quantitative PCR or after one round of 
amplification (WGA2; Sigma-Aldrich) and applied to a general extended 
NimbleGen MM8 Deluxe Promoter HX1 array (MOgene, LC). 

Reporter assay 

Plasmids expressing pSatbl-luc construct bearing a promoter region con- 
taining the far-upstream CpG island and the p63 consensus sequence of 
mouse Satbl were generated. For reporter assay, HaCaT cells were 
seeded into 96-well plates (1 0 4 cells per well) 1 d before transfection. Cells 
were cotransfected with pANp63 (Antonini et al., 2006) or pcDNA3 plas- 
mids (120 ng/well; Invitrogen) along with the reporter pSatbl-Luc con- 
struct using Lipofectamine 2000 (0.5 pl/well; Invitrogen). After 48 h, 
reporter activity was detected using the Dual-Glo Luciferase Assay System 
(Promega). Two independent assays were performed in triplicate, and re- 
sults were normalized to the pNull-Renilla construct activity, the data were 
pooled, a mean ± SEM was calculated, and statistical analysis was per- 
formed using a Student's ftest. 

Statistical analysis and gene ontology analysis 

For statistical evaluation of differences in radial positions, intergene dis- 
tances, and 3D FISH signal volumes, normality of data distributions was 
evaluated by plotting them on histograms. Equality of their variances was 
tested by the Levene's test. Significance of the differences for the pairwise 
comparisons was analyzed by a two-tailed ftest (a = 0.05). For more than 
two-sample comparisons, a one-way analysis of variance (ANOVA) fol- 
lowed by the Newman-Keuls test (a = 0.05) was performed. Enrichment 
of the Satbl target genes in the distinct genomic loci was assessed by 
using the hypergeometric or Fisher's exact tests. Functional annotation of 
the overrepresented and underrepresented genes was performed as previ- 
ously described (Sharov et al., 2009; Fessing et al., 2010) using the 
National Institute on Aging Array Analysis software (National Institutes of 
Health) according to previously published recommendations (Sharov 
etal., 2005a). 

Online supplemental material 

Fig. SI shows localization of the EDC in gene-dense regions of mouse 
chromosome 3. Fig. S2 shows validation of the effects of p63 knockdown 
after siRNA treatment in skin explants and keratinocytes as well as the data 
on the analyses of 3D FISH distances and gene expression in the epidermis 
of p63~/~, Satbl , and corresponding WT mice. Fig. S3 shows the 
results of the microarray, ChlP-on-chip, and bioinformatic analyses of the 
genes whose expression was altered in the skin of Satbl mice versus 
WT controls and which showed association of Satbl with the genomic 
regions within 200 kb from the transcription start sites in the ChlP-on-chip 
assay. Fig. S4 shows the results of ChlP-quantitative PCR analyses of 
keratinocytes demonstrating binding of Satbl to the promoter regions of 
the distinct genes within the EDC. Fig. S5 shows the results of 3D FISH 
analyses of radial positions of Rps27 and tor genes and differences in 
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the intergene distances in the EDC locus as well as the confocal stacks 
demonstrating intranuclear localization of the distinct EDC domains in the 
epidermis of Satbl and WT mice. Tables SI and S2 show the lists 
of the genes down- and up-regulated in the skin epithelium of the El 6.5 
p63~/~ versus WT mice, respectively. Table S3 shows the list of 3D FISH 
BAC probes generated to the distinct regions of mouse chromosome 3. 
Tables S4 and S5 show the list of genes that demonstrate Satbl binding to 
their regulatory regions and down- and up-regulated in the skin of the P0.5 
Satbl versus WT mice, respectively. Table S6 shows SATB1 binding 
sites in the EDC locus and its 5' and 3' flanking regions based on the ChlP- 
on-chip analysis. Table S7 shows the length of the distinct EDC domains in 
Satb 1 and WT mice. Tables S8 and S9 show the list of the primers used 
for quantitative RT-PCR and ChlP-quantitative PCR analyses, respectively. 
Table SI 0 shows the list of primary antibodies used for immunohistochemi- 
cal analyses. Online supplemental material is available at http://www 
.icb.org/cgi/content/full/jcb.201 101 148/DC1 . 
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